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Abstract. We give concisely the formulae governing diffusion of chemical ele- 
ments and their isotopes in quiescent stellar atmospheres, due to electrostatic, 
gravitational and radiation fields and to impacts between particles. Isotope 
segregation of heavy elements due to light-induced drift is emphasized. 
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1. Introduction 

The diffusive separation of chemical elements and their isotopes in stellar at- 
mospheres can occur only in the case of lacking macroscopic motions, i.e. if 
the stellar wind and the meridional circulation are extremely weak and there 
is no convective turbulence. These conditions hold only for CP stars, and over- 
abundances of heavy elements in their atmospheres can reach several orders of 
magnitude. 

2. Diffusion equations of plasma components 

The element diffusion is described with sufficient accuracy if we take into account 
the collisions between different atomic particles and the presence of external 
gravitation, radiation and electrostatic fields (Sapar & Aret, 1995). Momentum 
transfer from particles of types i to particles j is described by 

^^^^ + VP, = p,a, + ^^,P^V^ - UjPjV, . (1) 

i 

Here the collision frequency of particle j with particles i and total collision 
frequency for particles j are to be found from expressions 

<^ji{vjt)vjifjfidvidvj , i/j^^i/j^^ (2) 

i 

where rii is the number density, Vji the relative particle velocity and fi the 
undisturbed velocity distribution. Indices indicate species of particles — all ions 
are treated separately and electrons are included as one of them. 

In the time-independent approximation we obtain from Eq. (n]) a system of 
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algebraic linear equations for diffusion velocities Vj: 

^jPjVj =Fj+^ i^ijPiVi , where Fj = - VP,- + pjUj . (3) 

i 

Using a buffer-gas approximation we neglect interactions between heavy ele- 
ments as small admixtures and take into account for each of them separately 
only interaction with buffer gas particles. Calculations showed, that adequate 
buffer gas mixture for modelling consists of H, He, C, N and O. The diffusion 
equations for a particular element can be derived from the continuity eqation 

^ + V(p,y,)=0, where Pe = ^Pj, V, = -J2pjV, (4) 

where summing is carried out over all particle species belonging to clement e. 
The electron density is to be found from the condition of gas electroneutrality. 

The external acceleration of plasma particles due to gravity, electrostatic 
field, radiative acceleration both in spectral continua and lines and due to light - 
induced drift (Atutov & Shalagin, 1988; Nasyrov & Shalagin, 1993) has the form 

«J = 9j + —E + aj,rad ■ (5) 

rrij 

Electrons, since they are much lighter than ions, tend to drift up in the 
atmosphere, generating an electrostatic field which blocks their escape. The 
electrostatic field has almost no influence on diffusion of heavy elements, but 
for light ions its lifting force can even prevail upon the opposite-directed gravity. 

The electrostatic field can be found from the main system of linear equations 
for diffusion velocities and the condition of no electric current ^ ■ riiZiVi = . 

In numerous studies (see Gonzalez et al., 1995; Michaud & Proffitt, 1993; 
Alecian et al., 1993; Ryabchikova, 1992 and references therein) the radiative 
expelling force for the majority of elements has been computed both for stellar 
atmospheres and for their envelopes. For line-rich metals the expelling force 
highly exceeds the gravity, thus being the dominant factor for the formation of 
observed metal overabundances in the CP stellar atmospheres.. 

Some essential discrepancies between the observations and theory have re- 
mained (Wahlgren et al., 1995). They can be partly removed when the light- 
induced drift is also taken into account (see also LeBlanc & Michaud, 1993). 

In free-free and bound-free electron transitions the momentum of the photon 
is transferred to the generated ions. In bound-bound transitions it is transferred 
to particles in the upper state u and the radiation fiux acts on them with a force: 

fui = ^ j iT'i (^ui{y) du , (Jui{y) = aliW{u^, a) , (6) 

where ni is the number density of particles in the lower state, Fi, is the monochro- 
matic fiux, aui {i') is the photon absorption cross-section in the transition I u 
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which can be written using the Voigt function W{u,a). The parameters of the 
Voigt function are u = {v — and a = r„i/(47rAi/£)). 

To find the radiative force acting on a particular ion, we need to summarize 
over aU transitions. The effect of hght-induced drift can be taken into account 
by substituting the Voigt function W{u, a) in the radiative-force expression for 
spectral lines by 

T,./ , „ ^dW(u,a) dW(u,a) dW(—u,a) 

The expression q = cV2MkT/h V is the ratio of the mean momentum of the 
ion species studied to the momentum of the absorbed light quantum, hvjc^ 
being about 10^. The light-induced drift is generated due to the lower mobility 
of atomic particles in the excited (upper) states resulting from larger impact 
cross-sections than in the ground (lower) state, yielding thus lower diffusion 
rates (particle mobilities) for excited states. 

The uncompensated drift rate due to the difference of cross-sections a is given 
by the expression D\ = l — ai/du = ^ajou- The spontaneous transitions which 
take place before collisions reduce the drift rate by factor Di = \ — Aj i^A+Vy) = 
v^l {^A -h Vy), where A is the probability (frequency) of spontaneous transitions. 
Thus, the total reduction factor is 2? = D\Di = Ai^/{A + i^u)- The effective 
value of 2Dq reduces to about 10^ in the upper atmosphere and reaches to about 
10^ in the deeper layers. 

The most effective case for light-induced drift is if in blends, the fluxes in the 
blue and red wing of a spectral line differ essentially, especially for isotope split- 
ting of resonance spectral lines of heavy elements where the isotopic wavelength 
difference of spectral lines is of the order of the thermal Doppler width. In this 
case only the heaviest isotopes dominate in the atmospheres as observed, say, 
by HST for Hg, Tl and Pt in the spectrum of x Lup (Wahlgren et al., 1995). 
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